How stable is the photon? 
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Yes, the photon. While a nonzero photon mass has been under experimental and theoretical study 
for years, the possible implication of a finite photon lifetime lacks discussion. The tight experimental 
upper bound of the photon mass restricts the kinematically allowed final states of photon decay to 
the lightest neutrino and/or particles beyond the Standard Model. We discuss the modifications of 
the well-measured cosmic microwave background spectrum of free streaming photons due to photon 
mass and lifetime and obtain model-independent constraints on both parameters — most importantly 
a lower direct bound of 3yrs on the photon lifetime, should the photon mass be at its conservative 
upper limit. 

PACS numbers: 14.70.Bh,12.20.-m,12.60.-i,98.80.-k 
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Classical electrodynamics, as encoded in Maxwell's 
equations, can be readily extended to allow for a nonzero 
photon mass; the resulting Proca equations [l| then de- 
scribe the behavior of a massive spin-1 field, and have 
since been used to set an impressive upper limit on the 
photon mass of m < 2 x 10~^^kg [2|, or 10~^^cV in 
the natural units used in this letter (ft = c = fc^ = 1). 
A nonzero photon mass is often dismissed on theoret- 
ical grounds, as the insertion of a mass term to the 
Lagrangian of quantum electrodynamics (QED) breaks 
gauge invariance and might therefore spoil renormaliz- 
ability, i.e. the consistency of the theory at quantum level. 
This is however not the case, as the Proca Lagrangian 
can be viewed as a gauge-fixed version of the Stiickelberg 
Lagrangian [3|, which restores gauge invariance. For an 
exhaustive review we refer to [J|. To the point: gauge 
bosons of abelian symmetries are permitted a mass by 
means of the Stiickelberg mechanism — retaining gauge 
invariance, unitarity, and renormalizability. 

The question of a photon mass in QED is then purely 
experimental, as there is no theoretical prejudice against 
a small m over m = 0.^ However, we already know 
that QED is just the low-energy approximation of the 
Glashow-Weinberg-Salam model of electroweak inter- 
actions, so our above motivation for the photon mass 
might be in danger. Fortunately, the electroweak gauge 
group SU{2)l X U{1)y still features an abelian factor — 
the hypcrcharge U{1)y — that can be used in a Stiick- 
elberg mechanism. The resulting mass for the hypcr- 
charge gauge boson eventually generates again a mas- 
sive photon [6|.^ A detailed discussion of this procedure 
and its implications can be found in Ref. [J]. Since the 
Stiickelberg mechanism only works for abelian groups, 
the grand unification of the Standard Model (SM) gauge 
group SU{3)c X SU{2)l x U{1)y into a simple nonabelian 



^ A small m is technically natural [3] , in that all radiative correc- 
tions are proportional to m. 

^ The same trick works, for example, in simple left-right symmetric 
models Q, where the hypercharge U{1)y itself results from the 
breakdown of 5(7(2)^ X U{1)b—l: A Stiickelberg mass of the 
B — L boson trickles down and makes the photon massive. 



group like S'C/(5), S'O(IO) or Er would necessarily re- 
sult in a truly massless photon [8|. Turning this around, 
the discovery of a massive photon would exclude a huge 
number of grand unified theories — and, obviously, be a 
spectacular finding in its own right. 

Let us now move on to the key point of this letter: 
If one can constrain the mass of a photon, one should 
also be able to constrain its lifetime. Massless photons 
in QED are stable purely due to kinematical reasons, 
there are no additional quantum numbers that forbid a 
decay; a photon with mass above MeV would certainly 
decay into an electron-positron pair. Recalling the tight 
upper bound on the photon mass though, there are not 
many possible final states — indeed, only one known par- 
ticle could be even lighter than the photon: the lightest 
neutrino vi. This is because current neutrino-oscillation 
experiments can only fix the two mass-squared differences 
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neutrinos, so the absolute mass scale is not known as 
of yet [9|. Kinematically, this opens up the possibility 
of a decay 7 -^ i/ii/i — should mi < m/2 hold.'^ This 
loop-suppressed process can be calculated in the SM (us- 
ing e.g. a seesaw mechanism to make neutrinos massive 
in a renormalizable way), and is of course ridiculously 
small [12| — being suppressed by the small photon mass, 
the heavy particles in the loop and maybe the smallest 
neutrino mass, depending on the operator that induces 
this decay. We also note that one of the side effects of a 
massive hypercharge boson — besides a massive photon — 
are tiny electric charge shifts of the known (chiral) el- 
ementary particles [J, |6[. The neutrino then picks up 
a charge Qi/ oc ern^/M^^, which gives rise to a corre- 
spondingly small tree- level decay rate 7 — )- uivi. Still, 
unmeasurable small SM rates never stopped anyone from 



^ The naive prototype model — augmenting the SM by only two 
right-handed neutrinos (SM-|-2i/fl) — is problematic, as the ini- 
tially massless ui will unavoidably pick up a finite mass at loop 
level llOl . which can be too largo for our purposes llll . Fine- 
tuned solutions aside, we can obtain a simple valid model by 
imposing a, B — L symmetry on the SM-|-2!^fl, resulting in two 
Dirac neutrinos and one exactly massless Weyl neutrino. 



looking for a signal, as it would be a perfect sign for new 
physics. 

Particles beyond the SM could not only increase the 
rate 7 -^ J^iz^i, but also serve as final states themselves, 
as some extensions of the SM feature additional (close to) 
massless states; examples include sterile neutrinos, hid- 
den photons, Goldstonc bosons and axions. These weakly 
interacting sub-cV particles 13| are less constrained than 
neutrinos, and photon decay might be an indirect effect 
of these states. Although mainly of academic interest, we 
also mention that a massive photon provides the possibil- 
ity of faster-than-light particles — and a decaying photon 
even predicts them. The question of photon decay is 
therefore obviously relevant even if the lightest neutrino 
turns out to be an inaccessibly heavy final state. 

Following the above motivation, we set out to find 
limits on the photon mass m and lifetime r^ as model- 
independent parameters. Most importantly, we do not 
care about the daughter particles for now. Due to the 
small allowed values for m, all measureable photons 
around us are highly relativistic, making a decay hard 
to observe because of time dilation. Correspondingly, 
a good limit on r^ needs a large number of low-energy 
photons from well-known far-away sources. Seeing as we 
have access to very accurate measurements of the cosmic 
microwave background (CMB) — consisting of the oldest 
photons in the visible Universe — we will take m and T-y 
as parameters that will modify the black-body radiation 
law — given by the Planck spectrum — and fit the CMB 
spectrum to obtain bounds on both parameters. Similar 
analyses have been performed to obtain a limit on the 
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neutrino lifetime in the channels f^ — >■ 7 i/j 
our case, we are however not looking for a spectral line on 
top of the CMB, but rather a diminished overall intensity 
and change of shape. 

Before delving into the details, let us present a back- 
of-the-envelope estimate: CMB photons with low ener- 
gies around meV have a lifetime r = 7LT-y that is in- 
creased by a relativistic Lorentz factor 7l — E/m ~ 
ImeV/lO-i^eV = lO^^. This lifetime has to be com- 
pared to the age of the Universe tg — 13.8 x lO^yrs, 
further taking into account the measured CMB accuracy 
of about 10^''. We therefore expect a lifetime constraint 
in the ballpark of years, which will be confirmed by the 
more refined analysis below. 

The photon mass changes the spectral energy density 
of black-body radiation to 
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because of the modified dispersion relation p^ ^ E^ — m^ , 
but it is unclear how to include the decay width. The 
expansion of the Universe also needs to be taken into 
account, as the black-body spectrum no longer stays in 
shape for m ^ 0. Let us therefore give a brief deriva- 
tion of the energy spectrum of massive unstable photons 
during cosmic expansion. 

Ignoring the width for a moment, the number density 



of massive photons right after decoupling (at the time of 
last scattering t^ ~ 400 000 years) is given by [16| 
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where p — p^ a{tL)/a{t) is the redshifted momentum, 
T the temperature at time i, and g the number of spin 
states. We take g = 2, because only the transverse modes 
are excited before decoupling (this implicitly constrains 
m, as discussed below). The chemical potential of mass- 
less photons is zero, and since we assume that as our 
initial condition at t^, we set it to zero in all our calcu- 
lations. 

Including the width, we can write down the differential 
equation for the time evolution of the number density 
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The first term on the right-hand side describes the num- 
ber density dilution due to the expansion of the Universe, 
while the second one is due to photon decay. The width 
can be obtained from the rest- frame width Fq = 1/x.y 
by a Lorentz boost: T{p) ~ Fo — . We use the bound- 
ary condition n(p, t^) — no{p, t^) and obtain the number 
density today 
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The integral can be evaluated to 
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with the comoving distance of the surface of last scatter- 
ing dL = J^° a{to)/a{t) dt ~ 47 billion lightyears. Overall 
we have: 
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The energy density relevant for the CMB spectrum 
is then obtained by multiplying n{p,to) with E = 



vV + ' 
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Figure 1: CMB spectral distribution for Fotom — (gray), 
Totom = 2 X 10~^ meV (dashed red) and Totom = 10"'' meV 
(dotted blue) using Eq. (O, as well as the COBE data (error- 
bars multiplied by 1000 to be visible). In all cases the mass 
is m < lO^*" eV and has no visible effect. 



where we approximated 
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because a(ti)/a(^o) — 8 x 10^*. Due to this approxima- 
tion, the limit p{E — > m, T) is nonzero, which is however 
of no importance for the CMB analysis in this letter. 

Eq. ([7]) is the key equation of this letter and will now 
be used to set constraints on m and Tq. For illustra- 
tive purposes we show the spectrum for various values 
in Fig. [TJ As expected from time-dilation arguments, the 
low energy part of the spectrum shows the strongest devi- 
ations, which fortunately also features the smallest error 
bars. 

Using the COBE (COsmic Background Explorer) 
dataset of the CMB [17[ we can construct a simple x^ 
function to fit the spectrum from Eq. ^.^ The best 
fit values are at m = = Fq, so we can only obtain 
exclusion ranges, shown in Fig. [51 The limit on the pho- 
ton mass is not competitive with other experiments — 
TO < 3 X 10~^ eV — but for the photon width we find the 
only existing (and model-independent) bound 
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at 95% C.L. This would correspond to a photon lifetime 
of only three years, should the photon mass be close to 
its current bound! Another useful form of the constraint 
is given by 
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Figure 2: Constraints on photon mass m and lifetime r-y from 
the CMB spectrum. 



For two-particle fermionic final states X, the decay rate 
7 — )► XX from (effective) interactions like gXj^^XA^ will 
be of the form Fq - g'^m/4Tr ^. With Eq. ^ we 
can constrain g < 0.03 e, which corresponds to a very 
large effective electric charge and is excluded by other 
experiments [19[.^ In particular, final state neutrinos are 
far better constrained by their electric properties (see 
e.g. Ref. [15[ for a recent review) to be relevant in pho- 
ton decay. Our complementary and model-independent 
approach should be interesting nonetheless, as it consti- 
tutes the only direct constraint on the photon lifetime as 
of yet. 

Let us make a couple more comments to illustrate that 
our analysis above is somewhat inconsistent. Our ap- 
proach basically assumed a vanishing/negligible number 
density of Stiickelberg scalars (j) and daughter particles 
X prior to photon decoupling. To ensure this, m and Fq 
need to be small: has only the interaction mA'^d^cl), so 
for small mass m, it will not be in equilibrium with the 
rest of the SM. The creation rate of (j) via ej ^^ e(j) is pro- 
portional to c?vr? jT ^ which has to be smaller than the 
expansion rate of the Universe H{T) ~ T^/Mpi — at least 
before weak decoupling around T ~ 1 MeV — in order to 
not put (j) in thermal equilibrium at Big Bang nucleosyn- 
thesis (BBN). For m < 10"'^ eV, only the transversal po- 
larizations of the photon are excited, making it okay to 
treat the photon as massless before BBN. For the initial 
condition of our black-body calculation however, we need 
to ensure that only the two transverse degrees of freedom 
of the photon are excited at the surface of last scattering 
at T ~ 0.25 eV. This requires m < 5 x 10~i^eV, mak- 
ing our approach a little inconsistent, because at these 
low masses the primordial plasma — consisting mainly of 
partly ionized hydrogen and helium — can not be ignored. 
We will remark on this below. 

On to the daughter particles: the interaction rate of 
photons with their will-be daughter particles at tempera- 
ture T will be something like ToT/m, as it should be finite 
in the limit m — > 0. This rate has to be smaller than the 



"^ Ground based and balloon experiments probe the CMB down 
to energies ~ 10~® eV, which typically have much larger errors. 
Additionally, there is an excess at low energies that is not un- 
derstood yet [\M . so wc do not include that data. 



^ It is of course trivial to reinterpret bounds on milli-chargcd par- 
ticles [l9| in terms of photon decay. 



expansion rate of the Universe at BBN — unless the final 
daughter particles are neutrinos. This gives the condition 
Fq -C IQ^^^m < 10^'**' eV, which is far stronger than the 
bound we obtained from the CMB analysis above, di- 
rectly related to the fact that the mini-charge of new ul- 
tralight particles is tightly constrained [l9|. One should 
be careful with the above constraint though, because ad- 
ditional degrees of freedom at BBN are still allowed by 
cosmological observations [20[ ■ 

Having discussed the initial conditions of our 
analysis — which degrees of freedom are present at 
recombination — it is time to scrutinize our main assump- 
tion: that the photons are free streaming. This is usually 
a very good approximation, as the density of ionized hy- 
drogen is rather small after recombination, but it is still 
large enough to induce a plasma mass as large as 10~^ eV 
to the photon. Further complications arise from the non- 
ionized hydrogen and helium, as they effectively make 
the Universe a refractive medium — changing the disper- 
sion relation of on-shell pho tons even further. This has 
been emphasized in Ref. [21| , where CMB constraints on 
photon oscillations into hidden photons [22| have been 
discussed. Their analysis (and phenomenology) is very 
similar to our discussion of photon decay, but in our 
case the inclusion of the plasma is more difficult. The 
photon in a medium requires a careful treatment, as it 
becomes just one of several quasiparticles that can be ex- 
cited. (A well-studied example relevant to our discussion 
is the decay of plasmons — effectively massive photons — 
into neutrinos as a mechanism to cool stars 15|,|23|.) This 
makes it difficult, if not impossible, to constrain the prop- 
erties of a free photon — namely m and r^ — through a 
study of these quasiparticles; certainly not in the model- 
independent way we aspired to. Naively reinterpreting r^ 
as an effective coupling of the daughter particles to the 
photons — and further ignoring the vacuum mass m in the 
dense plasma — would lead back to the usual bounds on 
milli-charged particles [19[ . 

In conclusion, a massive photon sounds crazy and ex- 
otic, but it really is not. A massless photon is neither a 
theoretical prediction nor a necessity, but rather a phe- 
nomenological curiosity. We should try to understand 
why this parameter in the Lagrangian (that we can just 
write down) is so small. This is similar to the strong 
CP problem [2J|, and in both cases experiments so far 
have only come up with upper bounds for these param- 
eters. Independent of its actual value, a nonzero pho- 
ton mass immediately opens up the possibility of photon 
decay — even in the SM — which can, and should, also be 
constrained. Using the long-lived low-energy photons of 
the cosmic microwave background, we were able to derive 
the first direct bound on the photon lifetime in this let- 
ter. Using the largest allowed value for the photon mass 
from other experiments, we find a lower limit of about 
3yrs on the photon lifetime. 

A proper study of the challenging, but important, ef- 
fects of the primordial plasma on this limit lies outside 
the scope of this letter and will be left for future work. 
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